 JOINTMAGINV.EXE Version 5.00.04
 -------------------------------

 *** Ridge Regression Iteration    1.
 Return status code irr_lsstat =           0
 Optimum damping factor =   47.3500138356776     
 Optimum ridge parameter beta =   2242.02381023886     
 Damping factor range qa,qb   =  -133.932864491883        133.932895395110     
 kbest, beta_best =           5   2242.02381023886     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.012782466314       0.000870680       0.000796000       0.000074680       9.3818552  0.13393290E+02  0.20691366E-01          0          0
     2     _Null_     -0.012743254978     -19.250473439       0.000000000     -19.250473439      -5.3473537  0.42792788E+01  0.21080419E-01          0          0
     3     Ztcd       -0.011356587677     304.124109337     318.678636000     -14.554526663      -3.6506914  0.31215381E-05  0.20747651E-01          0          0
     4     Rad_a       0.012782488043     160.406639703     155.681248000       4.725391703       3.0352992  0.15190679E+01  0.20691364E-01          0          0
 Initial best average relative RMS error (%) =   109.832805727054     
 Current best average relative RMS error (%) =   99.6237051805121     
 qmax =   13.3932895707264     
 qmed =   2.89917335781978     
 qmin =  3.121538055503796E-006
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947048587786  0.17116417D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319509168712    0.302590707013  0.66197099D-03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000233067    0.000000220726  0.78027874D-03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.113420074325    0.107414321216  0.27050642D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  1
 qa,qb   =  -133.932864491883        133.932895395110     
 
 Summary Table - Ridge Regression Iteration    1
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.93289571     0.00332143     0.00166072     0.06019215     0.00000000     102.21676939         T & T         T & T
     2   108.22860265     0.00506000     0.00253000     0.06009399     0.02260757     102.02310991         T & T         T & T
     3    85.23002465     0.00808502     0.00404251     0.05992163     0.06245794     101.67974410         T & T         T & T
     4    64.93716172     0.01369084     0.00684542     0.05959688     0.13808349     101.02162407         T & T         T & T
     5    47.35001384     0.02486254     0.01243127     0.05892847     0.29606624      99.62370518         T & T         T & T
     6    32.46858101     0.04883177     0.02441589     0.05739378     0.67174482      96.22658557         T & F         F & F
     7    20.29286325     0.10209341     0.05104671     0.05345273     1.74296866      86.98016348         T & F         F & F
     8    10.82286055     0.20563375     0.10281687     0.04502167     5.67883655      79.94151499         T & F         F & F
     9     4.05857290     0.34807484     0.17403742     0.17646376    23.84721088      78.46110351         T & F         F & F
    10     3.00000000     0.39690944     0.19845472     0.34993115    32.34177973      67.71834566         T & F         F & F
    11     1.00000000     0.70664022     0.35332011     1.07689881    58.63937658      82.92459828         T & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   9.29512860839830     
 Total reduction in av_qd_rel_rmse (%) =   9.29512860839830     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   2242.02381023886     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.85638997E-05  0.29264141E-02  0.87067957E-03  0.17097119E-04       0.0196365  0.10000000E+01        1            0
     2         _Null_     0.88966533E-05  0.29827258E-02 -0.19250473E+02  0.45058256E+01       0.2340631  0.31950917E+00        1            0
     3         Ztcd       0.86121576E-05  0.29346478E-02  0.30412411E+03  0.37610249E+01       0.0123667  0.23306732E-06        1            0
     4         Rad_a      0.85638980E-05  0.29264138E-02  0.16040664E+03  0.10818279E+01       0.0067443  0.11342007E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.29264141E-02  0.87067957E-03  0.17097119E-04       1.0000000     171.1641711     0.020691366     0.001062766     2.482328905
     2         _Null_     0.29827258E-02 -0.19250473E+02  0.45058256E+01       0.3195092       0.0006620     0.021080419     0.000689526     1.551644740
     3         Ztcd       0.29346478E-02  0.30412411E+03  0.37610249E+01       0.0000002       0.0007803     0.020747651     0.001008523     2.342869440
     4         Rad_a      0.29264138E-02  0.16040664E+03  0.10818279E+01       0.1134201       0.0027051     0.020691364     0.001062768     2.482333822
 Total relative Backus-Gilbert R(spread) =   8.85917690655341     

 *** Ridge Regression Iteration    2.
 Return status code irr_lsstat =           0
 Optimum damping factor =   47.3515685725267     
 Optimum ridge parameter beta =   2242.17104627870     
 Damping factor range qa,qb   =  -133.937251343266        133.937292778428     
 kbest, beta_best =           5   2242.17104627870     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.014531166271       0.000941394       0.000870680       0.000070715       8.1218022  0.13393729E+02  0.20689940E-01          0          0
     2     _Null_     -0.015283007630     -36.193124677     -19.250473439     -16.942651238      -4.7062920  0.42798564E+01  0.21079175E-01          0          0
     3     Ztcd       -0.012295787178     292.505103680     304.124109337     -11.619005657      -3.0248051  0.41853699E-05  0.20745876E-01          0          0
     4     Rad_a       0.014531212603     164.621538620     160.406639703       4.214898917       2.6276337  0.15135534E+01  0.20689938E-01          0          0
 Initial best average relative RMS error (%) =   99.6237051805121     
 Current best average relative RMS error (%) =   96.4978906196726     
 qmax =   13.3937293196965     
 qmed =   2.89670489966486     
 qmin =  4.185369913719807E-006
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947079682733  0.20548954D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319541804707    0.302631551022  0.90204345D-03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000312487    0.000000295950  0.10582478D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.113004625121    0.107024384507  0.34475827D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  2
 qa,qb   =  -133.937251343266        133.937292778428     
 
 Summary Table - Ridge Regression Iteration    2
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.93729320     0.00379359     0.00189680     0.05371991     1.59958533      99.22422877         T & T         T & T
     2   108.23215620     0.00577950     0.00288975     0.05360313     1.63114150      99.01665706         T & T         T & T
     3    85.23282310     0.00923516     0.00461758     0.05339825     1.68676238      98.65026371         T & T         T & T
     4    64.93929389     0.01564012     0.00782006     0.05301289     1.79230471      97.95371550         T & T         T & T
     5    47.35156857     0.02840894     0.01420447     0.05222244     2.01273227      96.49789062         T & T         T & T
     6    32.46964715     0.05582912     0.02791456     0.05042288     2.53659362      93.08276763         T & F         F & F
     7    20.29352963     0.11693962     0.05846981     0.04592721     4.02754387      84.54078348         T & F         F & F
     8    10.82321600     0.23779415     0.11889708     0.03818923     9.45657058      77.07010477         T & F         F & F
     9     4.05870626     0.43946938     0.21973469     0.22729944    32.70779219      68.69919817         T & F         F & F
    10     3.00000000     0.52948079     0.26474039     0.43821400    42.24036122      69.53560416         T & F         F & F
    11     1.00000000     1.10393641     0.55196820     1.16128706    66.49132935      88.00329168         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   3.13762126712275     
 Total reduction in av_qd_rel_rmse (%) =   12.1411039434975     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   2242.17104627870     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.85627178E-05  0.29262122E-02  0.94139444E-03  0.14240200E-04       0.0151267  0.10000000E+01        1            0
     2         _Null_     0.88956051E-05  0.29825501E-02 -0.36193125E+02  0.33064373E+01       0.0913554  0.31954180E+00        1            0
     3         Ztcd       0.86106730E-05  0.29343948E-02  0.29250510E+03  0.27728806E+01       0.0094798  0.31248727E-06        1            0
     4         Rad_a      0.85627156E-05  0.29262118E-02  0.16462154E+03  0.84877204E+00       0.0051559  0.11300463E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.29262122E-02  0.94139444E-03  0.14240200E-04       1.0000000     205.4895366     0.020689940     0.001062789     2.482724047
     2         _Null_     0.29825501E-02 -0.36193125E+02  0.33064373E+01       0.3195418       0.0009020     0.021079175     0.000689376     1.551488882
     3         Ztcd       0.29343948E-02  0.29250510E+03  0.27728806E+01       0.0000003       0.0010582     0.020745876     0.001008883     2.344106031
     4         Rad_a      0.29262118E-02  0.16462154E+03  0.84877204E+00       0.1130046       0.0034476     0.020689938     0.001062791     2.482730624
 Total relative Backus-Gilbert R(spread) =   8.86104958359634     

 *** Ridge Regression Iteration    3.
 Return status code irr_lsstat =           0
 Optimum damping factor =   47.3529282115652     
 Optimum ridge parameter beta =   2242.29981020965     
 Damping factor range qa,qb   =  -133.941085086879        133.941138278662     
 kbest, beta_best =           5   2242.29981020965     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.013513112495       0.000997730       0.000941394       0.000056336       5.9842670  0.13394114E+02  0.20688694E-01          0          0
     2     _Null_     -0.014564509016     -48.601506837     -36.193124677     -12.408382160      -3.4467728  0.42803305E+01  0.21078081E-01          0          0
     3     Ztcd       -0.011212445877     284.329940840     292.505103680      -8.175162840      -2.1946445  0.53729074E-05  0.20744332E-01          0          0
     4     Rad_a       0.013513173330     167.808739780     164.621538620       3.187201160       1.9360779  0.15088024E+01  0.20688691E-01          0          0
 Initial best average relative RMS error (%) =   96.4978906196726     
 Current best average relative RMS error (%) =   93.8059010050240     
 qmax =   13.3941138815953     
 qmed =   2.89456640236002     
 qmin =  5.372907400318411E-006
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947106875366  0.23986827D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319568020061    0.302665068947  0.11737637D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000401139    0.000000379922  0.13715257D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112646671847    0.106688437394  0.42398244D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  3
 qa,qb   =  -133.941085086879        133.941138278662     
 
 Summary Table - Ridge Regression Iteration    3
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94113882     0.00353990     0.00176995     0.04884893     0.41494138      96.14274043         T & T         T & T
     2   108.23526379     0.00539309     0.00269655     0.04871396     0.45009413      95.96078224         T & T         T & T
     3    85.23527035     0.00861798     0.00430899     0.04847717     0.51213388      95.64133067         T & T         T & T
     4    64.94115849     0.01459573     0.00729787     0.04803169     0.63003560      95.03988835         T & T         T & T
     5    47.35292821     0.02651516     0.01325758     0.04711760     0.87663947      93.80590101         T & T         T & T
     6    32.47057951     0.05212354     0.02606177     0.04503458     1.46321573      91.01552330         T & F         F & F
     7    20.29411240     0.10928627     0.05464314     0.03981453     3.13160787      84.43633173         T & F         F & F
     8    10.82352686     0.22334285     0.11167143     0.03054622     9.17320647      75.47713219         T & F         F & F
     9     4.05882291     0.42899772     0.21449886     0.23163952    34.04694798      70.78132625         T & F         F & F
    10     3.00000000     0.52663476     0.26331738     0.45006651    43.69343719      73.14192491         T & F         F & F
    11     1.00000000     1.14257461     0.57128731     1.17101842    67.10763297      85.80281469         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   2.78968752307608     
 Total reduction in av_qd_rel_rmse (%) =   14.5920926046982     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   2242.29981020965     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.85616846E-05  0.29260357E-02  0.99773000E-03  0.12198511E-04       0.0122263  0.10000000E+01        1            0
     2         _Null_     0.88946826E-05  0.29823954E-02 -0.48601507E+02  0.25408823E+01       0.0522799  0.31956802E+00        1            0
     3         Ztcd       0.86093819E-05  0.29341748E-02  0.28432994E+03  0.21393510E+01       0.0075242  0.40113944E-06        1            0
     4         Rad_a      0.85616817E-05  0.29260352E-02  0.16780874E+03  0.69013121E+00       0.0041126  0.11264667E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.29260357E-02  0.99773000E-03  0.12198511E-04       1.0000000     239.8682702     0.020688694     0.001062808     2.483069048
     2         _Null_     0.29823954E-02 -0.48601507E+02  0.25408823E+01       0.3195680       0.0011738     0.021078081     0.000689251     1.551368439
     3         Ztcd       0.29341748E-02  0.28432994E+03  0.21393510E+01       0.0000004       0.0013715     0.020744332     0.001009190     2.345167369
     4         Rad_a      0.29260352E-02  0.16780874E+03  0.69013121E+00       0.1126467       0.0042398     0.020688691     0.001062812     2.483077477
 Total relative Backus-Gilbert R(spread) =   8.86268233324706     

 *** Ridge Regression Iteration    4.
 Return status code irr_lsstat =           0
 Optimum damping factor =   32.4712848459861     
 Optimum ridge parameter beta =   1054.38433954917     
 Damping factor range qa,qb   =  -133.943983429737        133.944047327295     
 kbest, beta_best =           6   1054.38433954917     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.024782438105       0.001089999       0.000997730       0.000092269       9.2479280  0.13394405E+02  0.40699618E-01          0          0
     2     _Null_     -0.027820318489     -68.138484142     -48.601506837     -19.536977305      -5.4269381  0.42806685E+01  0.42310460E-01          0          0
     3     Ztcd       -0.019840816375     272.370540964     284.329940840     -11.959399876      -3.2825740  0.64542988E-05  0.40930020E-01          0          0
     4     Rad_a       0.024782592497     172.829512928     167.808739780       5.020773148       2.9919617  0.15052564E+01  0.40699601E-01          0          0
 Initial best average relative RMS error (%) =   93.8059010050240     
 Current best average relative RMS error (%) =   89.0857041317492     
 qmax =   13.3944047972724     
 qmed =   2.89296249414665     
 qmin =  6.454298798789365E-006
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947127446211  0.26858797D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319586320199    0.302688975294  0.14239827D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000481865    0.000000456388  0.16590144D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112379495070    0.106437704172  0.49360112D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  4
 qa,qb   =  -133.943983429737        133.944047327295     
 
 Summary Table - Ridge Regression Iteration    4
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94404797     0.00332234     0.00166117     0.04332407     0.35297705      93.49127310         T & T         T & T
     2   108.23761465     0.00506173     0.00253087     0.04316768     0.39143839      93.33021154         T & T         T & T
     3    85.23712167     0.00808875     0.00404437     0.04289325     0.45950175      93.04838140         T & T         T & T
     4    64.94256905     0.01370027     0.00685013     0.04237679     0.58921323      92.52091862         T & T         T & T
     5    47.35395677     0.02489177     0.01244588     0.04131632     0.86112770      91.45095011         T & T         T & T
     6    32.47128485     0.04894869     0.02447434     0.03889574     1.50856508      89.08570413         T & T         T & T
     7    20.29455327     0.10274023     0.05137012     0.03279665     3.34845402      83.72445468         T & F         F & F
     8    10.82376205     0.21108888     0.10554444     0.02141774     9.96886233      75.87303076         T & F         F & F
     9     4.05891117     0.42081463     0.21040732     0.22998673    36.10457725      72.77816440         T & F         F & F
    10     3.00000000     0.52469489     0.26234745     0.45326791    45.68918243      75.02209895         T & F         F & F
    11     1.00000000     1.17115679     0.58557840     1.17360240    68.08987980      84.21669794         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   5.03187627079240     
 Total reduction in av_qd_rel_rmse (%) =   18.8897128303027     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   1054.38433954917     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.33150540E-04  0.57576505E-02  0.10899993E-02  0.21436740E-04       0.0196667  0.10000000E+01        1            0
     2         _Null_     0.35957049E-04  0.59964197E-02 -0.68138484E+02  0.42110199E+01       0.0618009  0.31958632E+00        1            0
     3         Ztcd       0.33554126E-04  0.57925923E-02  0.27237054E+03  0.34915866E+01       0.0128193  0.48186529E-06        1            0
     4         Rad_a      0.33150510E-04  0.57576480E-02  0.17282951E+03  0.11664576E+01       0.0067492  0.11237950E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.57576505E-02  0.10899993E-02  0.21436740E-04       1.0000000     268.5879720     0.040699618     0.004089749     2.468970932
     2         _Null_     0.59964197E-02 -0.68138484E+02  0.42110199E+01       0.3195863       0.0014240     0.042310460     0.002607736     1.456693031
     3         Ztcd       0.57925923E-02  0.27237054E+03  0.34915866E+01       0.0000005       0.0016590     0.040930020     0.003875809     2.313547464
     4         Rad_a      0.57576480E-02  0.17282951E+03  0.11664576E+01       0.1123795       0.0049360     0.040699601     0.004089764     2.468982326
 Total relative Backus-Gilbert R(spread) =   8.70819375264824     

 *** Ridge Regression Iteration    5.
 Return status code irr_lsstat =           0
 Optimum damping factor =   32.4724009074499     
 Optimum ridge parameter beta =   1054.45682069415     
 Damping factor range qa,qb   =  -133.948565797824        133.948650225415     
 kbest, beta_best =           6   1054.45682069415     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.021379619111       0.001157018       0.001089999       0.000067019       6.1485246  0.13394865E+02  0.40696534E-01          0          0
     2     _Null_     -0.025381119555     -81.427785504     -68.138484142     -13.289301362      -3.6914726  0.42811650E+01  0.42308060E-01          0          0
     3     Ztcd       -0.016360429244     264.984749309     272.370540964      -7.385791655      -2.0960298  0.85280395E-05  0.40925527E-01          0          0
     4     Rad_a       0.021379828026     176.267466256     172.829512928       3.437953328       1.9892166  0.14997384E+01  0.40696512E-01          0          0
 Initial best average relative RMS error (%) =   89.0857041317492     
 Current best average relative RMS error (%) =   85.5556025477431     
 qmax =   13.3948651078219     
 qmed =   2.89045169610134     
 qmin =  8.528039468991383E-006
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947159995082  0.31900891D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319612397659    0.302724076995  0.19098912D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000636665    0.000000603023  0.22151220D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111963683599    0.106047522007  0.62187662D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  5
 qa,qb   =  -133.948565797824        133.948650225415     
 
 Summary Table - Ridge Regression Iteration    5
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94865108     0.00289615     0.00144807     0.02978214     0.63383624      88.81716527         T & T         T & T
     2   108.24133437     0.00441261     0.00220631     0.02957235     0.67909922      88.69330183         T & T         T & T
     3    85.24005100     0.00705202     0.00352601     0.02920410     0.75941946      88.47748772         T & T         T & T
     4    64.94480096     0.01194621     0.00597311     0.02851066     0.91312044      88.07667330         T & T         T & T
     5    47.35558427     0.02171222     0.01085611     0.02708493     1.23697435      87.27547445         T & T         T & T
     6    32.47240091     0.04273231     0.02136615     0.02382035     2.01132762      85.55560255         T & T         T & T
     7    20.29525089     0.08993454     0.04496727     0.01550984     4.21120746      81.85835679         T & F         F & F
     8    10.82413420     0.18720264     0.09360132    -0.00140805    12.00992449      76.37866805         T & F         F & F
     9     4.05905086     0.40336350     0.20168175     0.21635036    40.05095774      75.03076479         T & F         F & F
    10     3.00000000     0.51646759     0.25823379     0.44600747    49.23568873      76.73596501         T & F         F & F
    11     1.00000000     1.20011091     0.60005546     1.16888981    69.55584348      81.73778596         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   3.96259042728727     
 Total reduction in av_qd_rel_rmse (%) =   22.1037813052343     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   1054.45682069415     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.33145478E-04  0.57572110E-02  0.11570182E-02  0.18047179E-04       0.0155980  0.10000000E+01        1            0
     2         _Null_     0.35952987E-04  0.59960810E-02 -0.81427786E+02  0.31394883E+01       0.0385555  0.31961240E+00        1            0
     3         Ztcd       0.33546569E-04  0.57919400E-02  0.26498475E+03  0.26147274E+01       0.0098675  0.63666483E-06        1            0
     4         Rad_a      0.33145440E-04  0.57572076E-02  0.17626747E+03  0.92577972E+00       0.0052521  0.11196368E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.57572110E-02  0.11570182E-02  0.18047179E-04       1.0000000     319.0089083     0.040696534     0.004089850     2.469406283
     2         _Null_     0.59960810E-02 -0.81427786E+02  0.31394883E+01       0.3196124       0.0019099     0.042308060     0.002607216     1.456567922
     3         Ztcd       0.57919400E-02  0.26498475E+03  0.26147274E+01       0.0000006       0.0022151     0.040925527     0.003877220     2.314898156
     4         Rad_a      0.57572076E-02  0.17626747E+03  0.92577972E+00       0.1119637       0.0062188     0.040696512     0.004089871     2.469421311
 Total relative Backus-Gilbert R(spread) =   8.71029367243231     

 *** Ridge Regression Iteration    6.
 Return status code irr_lsstat =           0
 Optimum damping factor =   20.2957197799996     
 Optimum ridge parameter beta =   411.916241388268     
 Damping factor range qa,qb   =  -133.951642383836        133.951743787737     
 kbest, beta_best =           7   411.916241388268     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.039154893888       0.001266733       0.001157018       0.000109715       9.4825681  0.13395174E+02  0.85463128E-01          0          0
     2     _Null_     -0.052224479579    -103.976786220     -81.427785504     -22.549000717      -6.2636113  0.42814652E+01  0.93844022E-01          0          0
     3     Ztcd       -0.026351537785     255.146230407     264.984749309      -9.838518902      -2.8518707  0.10242818E-04  0.86667232E-01          0          0
     4     Rad_a       0.039155519009     181.675138416     176.267466256       5.407672160       3.0678788  0.14961137E+01  0.85462989E-01          0          0
 Initial best average relative RMS error (%) =   85.5556025477431     
 Current best average relative RMS error (%) =   80.1048696214947     
 qmax =   13.3951744812019     
 qmed =   2.88878944730473     
 qmin =  1.024281830950159E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947181871083  0.35687809D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319627432766    0.302745309817  0.23160441D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000000764665    0.000000724277  0.26784050D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111690494650    0.105791211705  0.72407346D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  6
 qa,qb   =  -133.951642383836        133.951743787737     
 
 Summary Table - Ridge Regression Iteration    6
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95174481     0.00258093     0.00129047     0.01499018     0.45717716      85.34724663         T & T         T & T
     2   108.24383439     0.00393257     0.00196628     0.01472266     0.50752358      85.24976793         T & T         T & T
     3    85.24201980     0.00628547     0.00314273     0.01425299     0.59830277      85.08033662         T & T         T & T
     4    64.94630104     0.01064972     0.00532486     0.01336829     0.77472117      84.76703467         T & T         T & T
     5    47.35667812     0.01936379     0.00968189     0.01154832     1.15057402      84.14597452         T & T         T & T
     6    32.47315103     0.03814918     0.01907459     0.00737511     2.05297774      82.83500267         T & T         T & T
     7    20.29571978     0.08054850     0.04027425    -0.00329745     4.60884318      80.10486962         T & T         T & T
     8    10.82438436     0.17020894     0.08510447    -0.02593996    13.50463135      76.14920008         T & F         F & F
     9     4.05914478     0.39482932     0.19741466     0.19815191    42.67041747      75.66985849         T & F         F & F
    10     3.00000000     0.51559170     0.25779585     0.43178673    51.34862970      77.06063719         T & F         F & F
    11     1.00000000     1.22747414     0.61373707     1.15968524    70.06026926      80.02231582         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   6.37098303785147     
 Total reduction in av_qd_rel_rmse (%) =   27.0665361854056     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   411.916241388268     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.14666440E-03  0.12110508E-01  0.12667333E-02  0.33934580E-04       0.0267890  0.10000000E+01        1            0
     2         _Null_     0.18029105E-03  0.13427250E-01 -0.10397679E+03  0.57974935E+01       0.0557576  0.31962743E+00        1            0
     3         Ztcd       0.15156825E-03  0.12311306E-01  0.25514623E+03  0.45965065E+01       0.0180152  0.76466479E-06        1            0
     4         Rad_a      0.14666384E-03  0.12110485E-01  0.18167514E+03  0.16725492E+01       0.0092063  0.11169049E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.12110508E-01  0.12667333E-02  0.33934580E-04       1.0000000     356.8780863     0.085463128     0.017745733     2.429608927
     2         _Null_     0.13427250E-01 -0.10397679E+03  0.57974935E+01       0.3196274       0.0023160     0.093844022     0.010772509     1.223217358
     3         Ztcd       0.12311306E-01  0.25514623E+03  0.45965065E+01       0.0000008       0.0026784     0.086667232     0.016722598     2.226352394
     4         Rad_a      0.12110485E-01  0.18167514E+03  0.16725492E+01       0.1116905       0.0072407     0.085462989     0.017745847     2.429632455
 Total relative Backus-Gilbert R(spread) =   8.30881113397383     

 *** Ridge Regression Iteration    7.
 Return status code irr_lsstat =           0
 Optimum damping factor =   20.2964058765444     
 Optimum ridge parameter beta =   411.944091505424     
 Damping factor range qa,qb   =  -133.956137217714        133.956269937225     
 kbest, beta_best =           7   411.944091505424     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.028849235281       0.001335397       0.001266733       0.000068664       5.4205752  0.13395627E+02  0.85455889E-01          0          0
     2     _Null_     -0.045106046630    -118.661058119    -103.976786220     -14.684271899      -4.0789644  0.42818369E+01  0.93839713E-01          0          0
     3     Ztcd       -0.016007667615     250.621956599     255.146230407      -4.524273807      -1.3499402  0.13406011E-04  0.86653252E-01          0          0
     4     Rad_a       0.028850082344     184.861212935     181.675138416       3.186074519       1.7537207  0.14909884E+01  0.85455707E-01          0          0
 Initial best average relative RMS error (%) =   80.1048696214947     
 Current best average relative RMS error (%) =   77.1355227116559     
 qmax =   13.3956271277826     
 qmed =   2.88641265365295     
 qmin =  1.340601123872559E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947213878030  0.42014940D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319644380377    0.302771593127  0.30717251D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001000775    0.000000947948  0.35381739D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111304111414    0.105428799013  0.90550557D-02
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  7
 qa,qb   =  -133.956137217714        133.956269937225     
 
 Summary Table - Ridge Regression Iteration    7
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95627128     0.00198269     0.00099134    -0.02585405     0.77512833      79.92240563         T & T         T & T
     2   108.24749220     0.00302159     0.00151079    -0.02627259     0.82782941      79.86756157         T & T         T & T
     3    85.24490039     0.00483105     0.00241552    -0.02700724     0.92806689      79.77254136         T & T         T & T
     4    64.94849586     0.00819065     0.00409533    -0.02839050     1.13576397      79.59784873         T & T         T & T
     5    47.35827859     0.01491280     0.00745640    -0.03123376     1.60674905      79.25535584         T & T         T & T
     6    32.47424860     0.02947863     0.01473931    -0.03774063     2.78561657      78.54819172         T & T         T & T
     7    20.29640588     0.06289202     0.03144601    -0.05428982     6.15082281      77.13552271         T & T         T & T
     8    10.82475043     0.13900721     0.06950361    -0.08923484    17.30223464      75.17601599         T & F         F & F
     9     4.05928225     0.37844895     0.18922448     0.14966139    46.72119543      75.72322747         T & F         F & F
    10     3.00000000     0.50917467     0.25458734     0.38910775    54.01740062      76.74130290         T & F         F & F
    11     1.00000000     1.24551111     0.62275555     1.13292813    69.78381598      77.33889558         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   3.70682447130786     
 Total reduction in av_qd_rel_rmse (%) =   29.7700516698574     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   411.944091505424     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.14663855E-03  0.12109440E-01  0.13353975E-02  0.28821748E-04       0.0215829  0.10000000E+01        1            0
     2         _Null_     0.18027474E-03  0.13426643E-01 -0.11866106E+03  0.43710432E+01       0.0368364  0.31964438E+00        1            0
     3         Ztcd       0.15151461E-03  0.12309127E-01  0.25062196E+03  0.34789491E+01       0.0138813  0.10007752E-05        1            0
     4         Rad_a      0.14663782E-03  0.12109410E-01  0.18486121E+03  0.13373093E+01       0.0072341  0.11130411E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.12109440E-01  0.13353975E-02  0.28821748E-04       1.0000000     420.1494039     0.085455889     0.017746297     2.430097822
     2         _Null_     0.13426643E-01 -0.11866106E+03  0.43710432E+01       0.3196444       0.0030717     0.093839713     0.010770706     1.223124904
     3         Ztcd       0.12309127E-01  0.25062196E+03  0.34789491E+01       0.0000010       0.0035382     0.086653252     0.016728917     2.227912368
     4         Rad_a      0.12109410E-01  0.18486121E+03  0.13373093E+01       0.1113041       0.0090551     0.085455707     0.017746446     2.430128564
 Total relative Backus-Gilbert R(spread) =   8.31126365736344     

 *** Ridge Regression Iteration    8.
 Return status code irr_lsstat =           0
 Optimum damping factor =   20.2967586254151     
 Optimum ridge parameter beta =   411.958410698363     
 Damping factor range qa,qb   =  -133.958443086829        133.958596687683     
 kbest, beta_best =           7   411.958410698363     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.022684476172       0.001384926       0.001335397       0.000049528       3.7088648  0.13395860E+02  0.85452190E-01          0          0
     2     _Null_     -0.041208041220    -130.179595645    -118.661058119     -11.518537526      -3.1995938  0.42819774E+01  0.93837272E-01          0          0
     3     Ztcd       -0.009691029306     248.265838715     250.621956599      -2.356117885      -0.7126320  0.15515238E-04  0.86646327E-01          0          0
     4     Rad_a       0.022685493426     187.079442969     184.861212935       2.218230033       1.1999435  0.14884923E+01  0.85451981E-01          0          0
 Initial best average relative RMS error (%) =   77.1355227116559     
 Current best average relative RMS error (%) =   75.3831574842660     
 qmax =   13.3958598239207     
 qmed =   2.88523484251849     
 qmin =  1.551523774639842E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947230332132  0.45801236D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319649311639    0.302781523629  0.35775411D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001158211    0.000001097093  0.41131343D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111115847215    0.105252300863  0.10226844D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  8
 qa,qb   =  -133.958443086829        133.958596687683     
 
 Summary Table - Ridge Regression Iteration    8
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95859824     0.00164111     0.00082056    -0.06742398     1.04639596      77.03580588         T & T         T & T
     2   108.24937261     0.00250164     0.00125082    -0.06797602     0.96809145      77.00234460         T & T         T & T
     3    85.24638126     0.00400143     0.00200072    -0.06894425     0.83737586      76.94446718         T & T         T & T
     4    64.94962419     0.00678964     0.00339482    -0.07076465     0.63227674      76.83837676         T & T         T & T
     5    47.35910139     0.01238332     0.00619166    -0.07449533     0.61624547      76.63155603         T & T         T & T
     6    32.47481287     0.02458193     0.01229096    -0.08297311     1.87918318      76.20930099         T & T         T & T
     7    20.29675863     0.05311531     0.02655765    -0.10409619     5.91048725      75.38315748         T & T         T & T
     8    10.82493866     0.12326275     0.06163138    -0.14511076    18.52836670      74.26425827         T & F         F & F
     9     4.05935297     0.37757094     0.18878547     0.11157623    46.47815416      75.16144929         T & F         F & F
    10     3.00000000     0.51534763     0.25767382     0.35434682    52.82911298      76.03076612         T & F         F & F
    11     1.00000000     1.27071784     0.63535892     1.11148836    67.54512779      75.66702282         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   2.27180054764201     
 Total reduction in av_qd_rel_rmse (%) =   31.3655360206303     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   411.958410698363     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.14662536E-03  0.12108896E-01  0.13849256E-02  0.26437923E-04       0.0190898  0.10000000E+01        1            0
     2         _Null_     0.18026547E-03  0.13426298E-01 -0.13017960E+03  0.37529402E+01       0.0288289  0.31964931E+00        1            0
     3         Ztcd       0.15148811E-03  0.12308050E-01  0.24826584E+03  0.29923775E+01       0.0120531  0.11582114E-05        0            0
     4         Rad_a      0.14662451E-03  0.12108861E-01  0.18707944E+03  0.11840271E+01       0.0063290  0.11111585E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.12108896E-01  0.13849256E-02  0.26437923E-04       1.0000000     458.0123597     0.085452190     0.017746565     2.430344897
     2         _Null_     0.13426298E-01 -0.13017960E+03  0.37529402E+01       0.3196493       0.0035775     0.093837272     0.010769960     1.223103886
     3         Ztcd       0.12308050E-01  0.24826584E+03  0.29923775E+01       0.0000012       0.0041131     0.086646327     0.016731942     2.228671360
     4         Rad_a      0.12108861E-01  0.18707944E+03  0.11840271E+01       0.1111158       0.0102268     0.085451981     0.017746738     2.430380445
 Total relative Backus-Gilbert R(spread) =   8.31250058774219     

 *** Ridge Regression Iteration    9.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8250476393100     
 Optimum ridge parameter beta =   117.181656393331     
 Damping factor range qa,qb   =  -133.959774741155        133.959943448094     
 kbest, beta_best =           8   117.181656393331     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.035164006416       0.001457681       0.001384926       0.000072756       5.2534120  0.13395995E+02  0.17446803E+00          0          0
     2     _Null_     -0.101370919914    -155.890784431    -130.179595645     -25.711188786      -7.1419969  0.42820260E+01  0.22739399E+00          0          0
     3     Ztcd        0.009089845836     250.273045603     248.265838715       2.007206888       0.6114577  0.17041105E-04  0.18245930E+00          0          0
     4     Rad_a       0.035167550834     190.259318845     187.079442969       3.179875877       1.6997463  0.14871397E+01  0.17446658E+00          0          0
 Initial best average relative RMS error (%) =   75.3831574842660     
 Current best average relative RMS error (%) =   73.6041136812525     
 qmax =   13.3959945152205     
 qmed =   2.88458283537969     
 qmin =  1.704110496269767E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947239856245  0.48331520D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319649726206    0.302784960700  0.39426773D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001272105    0.000001204988  0.45286043D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111013757883    0.105156656058  0.11059410D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) =  9
 qa,qb   =  -133.959774741155        133.959943448094     
 
 Summary Table - Ridge Regression Iteration    9
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95994515     0.00140236     0.00070118    -0.10994821     2.06972245      75.32330010         T & T         T & T
     2   108.25046106     0.00213833     0.00106916    -0.11060937     1.96544980      75.30223971         T & T         T & T
     3    85.24723844     0.00342208     0.00171104    -0.11176764     1.78316181      75.26584635         T & T         T & T
     4    64.95027732     0.00581235     0.00290617    -0.11394067     1.44408899      75.19924835         T & T         T & T
     5    47.35957767     0.01062296     0.00531148    -0.11837433     0.79800933      75.06982044         T & T         T & T
     6    32.47513951     0.02119355     0.01059677    -0.12834568     1.23270909      74.80709575         T & T         T & T
     7    20.29696283     0.04646594     0.02323297    -0.15246176     5.75475402      74.29717227         T & T         T & T
     8    10.82504764     0.11327821     0.05663910    -0.19376645    19.12610591      73.60411368         T & T         T & T
     9     4.05939393     0.37846451     0.18923225     0.08109666    44.86419551      74.64099243         T & F         F & F
    10     3.00000000     0.52044618     0.26022309     0.32588436    50.53954402      75.42405780         T & F         F & F
    11     1.00000000     1.28477524     0.64238762     1.09378495    64.94944439      74.52063108         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   2.36000170646180     
 Total reduction in av_qd_rel_rmse (%) =   32.9853105417643     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.181656393331     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63370183E-03  0.25173435E-01  0.14576814E-02  0.52084923E-04       0.0357313  0.10000000E+01        1            0
     2         _Null_     0.12007738E-02  0.34652183E-01 -0.15589078E+03  0.87889980E+01       0.0563792  0.31964973E+00        1            0
     3         Ztcd       0.72481634E-03  0.26922413E-01  0.25027305E+03  0.59449692E+01       0.0237539  0.12721045E-05        0            0
     4         Rad_a      0.63368627E-03  0.25173126E-01  0.19025932E+03  0.22761726E+01       0.0119635  0.11101376E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25173435E-01  0.14576814E-02  0.52084923E-04       1.0000000     483.3151997     0.174468031     0.069770707     2.292141377
     2         _Null_     0.34652183E-01 -0.15589078E+03  0.87889980E+01       0.3196497       0.0039427     0.227393992     0.034977304     0.676438558
     3         Ztcd       0.26922413E-01  0.25027305E+03  0.59449692E+01       0.0000013       0.0045286     0.182459301     0.064232725     1.929409144
     4         Rad_a      0.25173126E-01  0.19025932E+03  0.22761726E+01       0.1110138       0.0110594     0.174466579     0.069771586     2.292208400
 Total relative Backus-Gilbert R(spread) =   7.19019747974984     

 *** Ridge Regression Iteration   10.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8250179383980     
 Optimum ridge parameter beta =   117.181013366639     
 Damping factor range qa,qb   =  -133.959390776892        133.959574075366     
 kbest, beta_best =           8   117.181013366639     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.025772139561       0.001509149       0.001457681       0.000051468       3.5307856  0.13395958E+02  0.17447060E+00          0          0
     2     _Null_     -0.096238450914    -178.414138582    -155.890784431     -22.523354151      -6.2564873  0.42818033E+01  0.22738934E+00          0          0
     3     Ztcd        0.018164338919     253.970291293     250.273045603       3.697245690       1.1194512  0.18514997E-04  0.18247131E+00          0          0
     4     Rad_a       0.025776078405     192.432926316     190.259318845       2.173607470       1.1424447  0.14881132E+01  0.17446903E+00          0          0
 Initial best average relative RMS error (%) =   73.6041136812525     
 Current best average relative RMS error (%) =   72.5812294418122     
 qmax =   13.3959575926866     
 qmed =   2.88495824735192     
 qmin =  1.851499735935787E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947237245428  0.50074487D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319633981639    0.302769212312  0.42728294D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001382133    0.000001309208  0.49129380D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111086738010    0.105225495716  0.11858663D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 10
 qa,qb   =  -133.959390776892        133.959574075366     
 
 Summary Table - Ridge Regression Iteration   10
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95957593     0.00117262     0.00058631    -0.17138414    14.58384807      73.55588223         T & T         T & T
     2   108.25016272     0.00178886     0.00089443    -0.17215130    14.46058298      73.54418733         T & T         T & T
     3    85.24700354     0.00286520     0.00143260    -0.17349255    14.24416096      73.52398602         T & T         T & T
     4    64.95009837     0.00487417     0.00243709    -0.17599946    13.83645405      73.48704246         T & T         T & T
     5    47.35944723     0.00893772     0.00446886    -0.18107548    12.99742393      73.41531321         T & T         T & T
     6    32.47505011     0.01797109     0.00898554    -0.19229136    11.07073306      73.26977842         T & T         T & T
     7    20.29690701     0.04025973     0.02012986    -0.21810775     6.10316846      72.98507131         T & T         T & T
     8    10.82501794     0.10450068     0.05225034    -0.25383466     7.60138756      72.58122944         T & T         T & T
     9     4.05938288     0.37850303     0.18925151     0.03977541    30.21539464      73.72736955         T & F         F & F
    10     3.00000000     0.52239051     0.26119525     0.28506999    35.85030505      74.46220884         T & F         F & F
    11     1.00000000     1.27960279     0.63980139     1.06654802    53.42055298      73.22171769         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   1.38971069452722     
 Total reduction in av_qd_rel_rmse (%) =   33.9166208480696     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.181013366639     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63372868E-03  0.25173968E-01  0.15091490E-02  0.50273043E-04       0.0333122  0.10000000E+01        1            0
     2         _Null_     0.12007294E-02  0.34651543E-01 -0.17841414E+03  0.81097417E+01       0.0454546  0.31963398E+00        1            0
     3         Ztcd       0.72495060E-03  0.26924907E-01  0.25397029E+03  0.54804084E+01       0.0215789  0.13821332E-05        0            0
     4         Rad_a      0.63371177E-03  0.25173632E-01  0.19243293E+03  0.21228053E+01       0.0110314  0.11108674E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25173968E-01  0.15091490E-02  0.50273043E-04       1.0000000     500.7448688     0.174470603     0.069769643     2.292038844
     2         _Null_     0.34651543E-01 -0.17841414E+03  0.81097417E+01       0.3196340       0.0042728     0.227389340     0.034980737     0.676532626
     3         Ztcd       0.26924907E-01  0.25397029E+03  0.54804084E+01       0.0000014       0.0049129     0.182471311     0.064225086     1.928925745
     4         Rad_a      0.25173632E-01  0.19243293E+03  0.21228053E+01       0.1110867       0.0118587     0.174469026     0.069770598     2.292111663
 Total relative Backus-Gilbert R(spread) =   7.18960887757566     

 *** Ridge Regression Iteration   11.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8249182114537     
 Optimum ridge parameter beta =   117.178854284663     
 Damping factor range qa,qb   =  -133.958149372852        133.958339145197     
 kbest, beta_best =           8   117.178854284663     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.021962270912       0.001552733       0.001509149       0.000043584       2.8879890  0.13395834E+02  0.17447731E+00          0          0
     2     _Null_     -0.094817390986    -199.952709594    -178.414138582     -21.538571012      -5.9829364  0.42815312E+01  0.22738584E+00          0          0
     3     Ztcd        0.021872252527     258.283972211     253.970291293       4.313680919       1.2916361  0.19168924E-04  0.18249541E+00          0          0
     4     Rad_a       0.021966417600     194.231191562     192.432926316       1.798265247       0.9344894  0.14900066E+01  0.17447568E+00          0          0
 Initial best average relative RMS error (%) =   72.5812294418122     
 Current best average relative RMS error (%) =   71.8069280557040     
 qmax =   13.3958341062089     
 qmed =   2.88576886050316     
 qmin =  1.916892367710625E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947228513615  0.50390605D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319616615222    0.302749971364  0.44022136D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001430962    0.000001355448  0.50704382D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111229099591    0.105359374676  0.12215338D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 11
 qa,qb   =  -133.958149372852        133.958339145197     
 
 Summary Table - Ridge Regression Iteration   11
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95834106     0.00109150     0.00054575    -0.20564737    18.60266725      72.55636787         T & T         T & T
     2   108.24916486     0.00166556     0.00083278    -0.20644786    18.47303356      72.54771483         T & T         T & T
     3    85.24621774     0.00266899     0.00133449    -0.20784569    18.24558140      72.53276570         T & T         T & T
     4    64.94949968     0.00454451     0.00227225    -0.21045248    17.81761690      72.50541791         T & T         T & T
     5    47.35901071     0.00834887     0.00417443    -0.21570675    16.93901510      72.45227031         T & T         T & T
     6    32.47475080     0.01686050     0.00843025    -0.22719578    14.93180182      72.34408392         T & T         T & T
     7    20.29671997     0.03820934     0.01910467    -0.25289306     9.80447361      72.12927427         T & T         T & T
     8    10.82491821     0.10214498     0.05107249    -0.28407122     3.47477175      71.80692806         T & T         T & T
     9     4.05934553     0.38102057     0.19051028     0.01580925    24.81821836      72.97848141         T & F         F & F
    10     3.00000000     0.52578412     0.26289206     0.26034179    30.66151687      73.71024020         T & F         F & F
    11     1.00000000     1.27758386     0.63879193     1.04889903    49.97512500      72.40654240         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =   1.06680665519577     
 Total reduction in av_qd_rel_rmse (%) =   34.6216027348407     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.178854284663     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63379419E-03  0.25175269E-01  0.15527331E-02  0.49960244E-04       0.0321757  0.10000000E+01        1            0
     2         _Null_     0.12006976E-02  0.34651084E-01 -0.19995271E+03  0.78712863E+01       0.0393657  0.31961662E+00        1            0
     3         Ztcd       0.72521532E-03  0.26929822E-01  0.25828397E+03  0.53111429E+01       0.0205632  0.14309616E-05        0            0
     4         Rad_a      0.63377667E-03  0.25174921E-01  0.19423119E+03  0.20609271E+01       0.0106107  0.11122910E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25175269E-01  0.15527331E-02  0.49960244E-04       1.0000000     503.9060504     0.174477315     0.069767704     2.291798838
     2         _Null_     0.34651084E-01 -0.19995271E+03  0.78712863E+01       0.3196166       0.0044022     0.227385843     0.034985152     0.676638835
     3         Ztcd       0.26929822E-01  0.25828397E+03  0.53111429E+01       0.0000014       0.0050704     0.182495409     0.064210935     1.927991452
     4         Rad_a      0.25174921E-01  0.19423119E+03  0.20609271E+01       0.1112291       0.0122153     0.174475681     0.069768693     2.291874252
 Total relative Backus-Gilbert R(spread) =   7.18830337696804     

 *** Ridge Regression Iteration   12.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8247949610836     
 Optimum ridge parameter beta =   117.176185949502     
 Damping factor range qa,qb   =  -133.956620456006        133.956813511476     
 kbest, beta_best =           8   117.176185949502     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.019840539329       0.001592271       0.001552733       0.000039538       2.5463194  0.13395682E+02  0.17448539E+00          0          0
     2     _Null_     -0.094081020861    -221.079978708    -199.952709594     -21.127269114      -5.8686859  0.42812479E+01  0.22738287E+00          0          0
     3     Ztcd        0.023700891346     262.896005085     258.283972211       4.612032873       1.3633613  0.19500553E-04  0.18252339E+00          0          0
     4     Rad_a       0.019844792070     195.831565413     194.231191562       1.600373851       0.8239531  0.14921904E+01  0.17448373E+00          0          0
 Initial best average relative RMS error (%) =   71.8069280557040     
 Current best average relative RMS error (%) =   71.1405241934220     
 qmax =   13.3956815461532     
 qmed =   2.88671916558725     
 qmin =  1.950055257547312E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947217725990  0.50181516D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319599110663    0.302729942831  0.44530611D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001455734    0.000001378897  0.51389251D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111393393248    0.105513796643  0.12400098D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 12
 qa,qb   =  -133.956620456006        133.956813511476     
 
 Summary Table - Ridge Regression Iteration   12
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95681546     0.00104986     0.00052493    -0.22795658    19.95943046      71.78894715         T & T         T & T
     2   108.24793206     0.00160228     0.00080114    -0.22877330    19.82714081      71.78167447         T & T         T & T
     3    85.24524691     0.00256834     0.00128417    -0.23019844    19.59513302      71.76910698         T & T         T & T
     4    64.94876002     0.00437555     0.00218777    -0.23285256    19.15895764      71.74610478         T & T         T & T
     5    47.35847137     0.00804762     0.00402381    -0.23818764    18.26497517      71.70135072         T & T         T & T
     6    32.47438098     0.01629476     0.00814738    -0.24978103    16.23004646      71.60994633         T & T         T & T
     7    20.29648885     0.03717676     0.01858838    -0.27528091    11.07732184      71.42621747         T & T         T & T
     8    10.82479496     0.10099720     0.05049860    -0.30391733     1.91902153      71.14052419         T & T         T & T
     9     4.05929933     0.38189047     0.19094523    -0.00277761    22.48994545      72.34696303         T & F         F & F
    10     3.00000000     0.52640922     0.26320461     0.24055271    28.41237603      73.09343747         T & F         F & F
    11     1.00000000     1.26896286     0.63448143     1.03404797    48.50592456      71.82277473         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.928049535505915     
 Total reduction in av_qd_rel_rmse (%) =   35.2283466469812     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.176185949502     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63387236E-03  0.25176822E-01  0.15922706E-02  0.50171505E-04       0.0315094  0.10000000E+01        1            0
     2         _Null_     0.12006713E-02  0.34650704E-01 -0.22107998E+03  0.77813223E+01       0.0351969  0.31959911E+00        1            0
     3         Ztcd       0.72552180E-03  0.26935512E-01  0.26289601E+03  0.52414681E+01       0.0199374  0.14557343E-05        0            0
     4         Rad_a      0.63385453E-03  0.25176468E-01  0.19583157E+03  0.20303443E+01       0.0103678  0.11139339E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25176822E-01  0.15922706E-02  0.50171505E-04       1.0000000     501.8151618     0.174485391     0.069765494     2.291514081
     2         _Null_     0.34650704E-01 -0.22107998E+03  0.77813223E+01       0.3195991       0.0044531     0.227382869     0.034989818     0.676746780
     3         Ztcd       0.26935512E-01  0.26289601E+03  0.52414681E+01       0.0000015       0.0051389     0.182523389     0.064194725     1.926913822
     4         Rad_a      0.25176468E-01  0.19583157E+03  0.20303443E+01       0.1113934       0.0124001     0.174483728     0.069766501     2.291590833
 Total relative Backus-Gilbert R(spread) =   7.18676551647979     

 *** Ridge Regression Iteration   13.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8246622168486     
 Optimum ridge parameter beta =   117.173312108869     
 Damping factor range qa,qb   =  -133.954975915912        133.955170598214     
 kbest, beta_best =           8   117.173312108869     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.018651811524       0.001629784       0.001592271       0.000037513       2.3559582  0.13395517E+02  0.17449400E+00          0          0
     2     _Null_     -0.093784119592    -242.090134888    -221.079978708     -21.010156180      -5.8361545  0.42809646E+01  0.22738022E+00          0          0
     3     Ztcd        0.024525102896     267.647479699     262.896005085       4.751474614       1.3856897  0.19664879E-04  0.18255279E+00          0          0
     4     Rad_a       0.018656119106     197.324521029     195.831565413       1.492955615       0.7623672  0.14944765E+01  0.17449232E+00          0          0
 Initial best average relative RMS error (%) =   71.1405241934220     
 Current best average relative RMS error (%) =   70.5270128832569     
 qmax =   13.3955172564701     
 qmed =   2.88772051524024     
 qmin =  1.966487894675500E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947206108955  0.49720622D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319581877163    0.302709906360  0.44637517D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001468019    0.000001390517  0.51615772D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111565417854    0.105675445340  0.12496098D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 13
 qa,qb   =  -133.954975915912        133.955170598214     
 
 Summary Table - Ridge Regression Iteration   13
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95517256     0.00102911     0.00051455    -0.24257783    20.69509070      71.12542375         T & T         T & T
     2   108.24660447     0.00157076     0.00078538    -0.24340748    20.56165556      71.11880448         T & T         T & T
     3    85.24420144     0.00251823     0.00125911    -0.24485464    20.32769832      71.10736444         T & T         T & T
     4    64.94796347     0.00429149     0.00214574    -0.24754792    19.88806327      71.08641965         T & T         T & T
     5    47.35789056     0.00789800     0.00394900    -0.25295419    18.98782784      71.04564101         T & T         T & T
     6    32.47398272     0.01601491     0.00800746    -0.26466568    16.94285987      70.96219683         T & T         T & T
     7    20.29623994     0.03667205     0.01833603    -0.29021865    11.79027706      70.79334143         T & T         T & T
     8    10.82466222     0.10046334     0.05023167    -0.31803727     1.03462982      70.52701288         T & T         T & T
     9     4.05924956     0.38205042     0.19102521    -0.01850093    21.18261775      71.77584299         T & F         F & F
    10     3.00000000     0.52584871     0.26292436     0.22339137    27.14799048      72.54620574         T & F         F & F
    11     1.00000000     1.25777283     0.62888641     1.02064024    47.68455148      71.35580703         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.862393575421239     
 Total reduction in av_qd_rel_rmse (%) =   35.7869332241918     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.173312108869     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63395537E-03  0.25178470E-01  0.16297839E-02  0.50639894E-04       0.0310715  0.10000000E+01        1            0
     2         _Null_     0.12006482E-02  0.34650371E-01 -0.24209013E+03  0.77626118E+01       0.0320650  0.31958188E+00        1            0
     3         Ztcd       0.72584349E-03  0.26941483E-01  0.26764748E+03  0.52196222E+01       0.0195019  0.14680194E-05        0            0
     4         Rad_a      0.63393737E-03  0.25178113E-01  0.19732452E+03  0.20148781E+01       0.0102110  0.11156542E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25178470E-01  0.16297839E-02  0.50639894E-04       1.0000000     497.2062238     0.174493997     0.069763192     2.291212464
     2         _Null_     0.34650371E-01 -0.24209013E+03  0.77626118E+01       0.3195819       0.0044638     0.227380216     0.034994522     0.676853556
     3         Ztcd       0.26941483E-01  0.26764748E+03  0.52196222E+01       0.0000015       0.0051616     0.182552795     0.064177786     1.925784823
     4         Rad_a      0.25178113E-01  0.19732452E+03  0.20148781E+01       0.1115654       0.0124961     0.174492318     0.069764208     2.291289901
 Total relative Backus-Gilbert R(spread) =   7.18514074477183     

 *** Ridge Regression Iteration   14.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8245269391915     
 Optimum ridge parameter beta =   117.170383457282     
 Damping factor range qa,qb   =  -133.953300971052        133.953496439001     
 kbest, beta_best =           8   117.170383457282     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.017931825262       0.001666270       0.001629784       0.000036487       2.2387371  0.13395350E+02  0.17450272E+00          0          0
     2     _Null_     -0.093700937890    -263.124994917    -242.090134888     -21.034860029      -5.8430167  0.42806859E+01  0.22737776E+00          0          0
     3     Ztcd        0.024838131359     272.459894125     267.647479699       4.812414426       1.3842800  0.19744237E-04  0.18258243E+00          0          0
     4     Rad_a       0.017936161106     198.754025093     197.324521029       1.429504064       0.7244432  0.14967736E+01  0.17450103E+00          0          0
 Initial best average relative RMS error (%) =   70.5270128832569     
 Current best average relative RMS error (%) =   69.9432080830336     
 qmax =   13.3953498413425     
 qmed =   2.88872973604560     
 qmin =  1.974423721454223E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947194270918  0.49146364D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319565065615    0.302690199336  0.44545549D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001473962    0.000001396128  0.51612619D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111738300144    0.105837877739  0.12547121D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 14
 qa,qb   =  -133.953300971052        133.953496439001     
 
 Summary Table - Ridge Regression Iteration   14
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95349841     0.00101843     0.00050921    -0.25300171    21.05501489      70.51321514         T & T         T & T
     2   108.24525162     0.00155454     0.00077727    -0.25384482    20.92105898      70.50692394         T & T         T & T
     3    85.24313607     0.00249244     0.00124622    -0.25531520    20.68622214      70.49605048         T & T         T & T
     4    64.94715176     0.00424827     0.00212413    -0.25805086    20.24504962      70.47614142         T & T         T & T
     5    47.35729870     0.00782117     0.00391058    -0.26353883    19.34213637      70.43737161         T & T         T & T
     6    32.47357687     0.01587165     0.00793583    -0.27541162    17.29343449      70.35798890         T & T         T & T
     7    20.29598628     0.03641593     0.01820797    -0.30123822    12.14572728      70.19698689         T & T         T & T
     8    10.82452694     0.10020007     0.05010004    -0.32928938     0.59135502      69.94320808         T & T         T & T
     9     4.05919884     0.38179529     0.19089764    -0.03287044    20.46831272      71.24004380         T & F         F & F
    10     3.00000000     0.52462658     0.26231329     0.20743825    26.43563302      72.03909409         T & F         F & F
    11     1.00000000     1.24519618     0.62259809     1.00780396    47.19988837      70.95256567         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.827774743827203     
 Total reduction in av_qd_rel_rmse (%) =   36.3184727731988     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.170383457282     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63403936E-03  0.25180138E-01  0.16662704E-02  0.51234997E-04       0.0307483  0.10000000E+01        1            0
     2         _Null_     0.12006272E-02  0.34650068E-01 -0.26312499E+03  0.77785702E+01       0.0295623  0.31956507E+00        1            0
     3         Ztcd       0.72616750E-03  0.26947495E-01  0.27245989E+03  0.52211059E+01       0.0191628  0.14739620E-05        0            0
     4         Rad_a      0.63402128E-03  0.25179779E-01  0.19875403E+03  0.20068172E+01       0.0100970  0.11173830E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25180138E-01  0.16662704E-02  0.51234997E-04       1.0000000     491.4636401     0.174502719     0.069760885     2.290907656
     2         _Null_     0.34650068E-01 -0.26312499E+03  0.77785702E+01       0.3195651       0.0044546     0.227377764     0.034999168     0.676958025
     3         Ztcd       0.26947495E-01  0.27245989E+03  0.52211059E+01       0.0000015       0.0051613     0.182582431     0.064160762     1.924649029
     4         Rad_a      0.25179779E-01  0.19875403E+03  0.20068172E+01       0.1117383       0.0125471     0.174501033     0.069761905     2.290985447
 Total relative Backus-Gilbert R(spread) =   7.18350015676209     

 *** Ridge Regression Iteration   15.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8243924508985     
 Optimum ridge parameter beta =   117.167471931067     
 Damping factor range qa,qb   =  -133.951636284866        133.951832102646     
 kbest, beta_best =           8   117.167471931067     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.017452612608       0.001702236       0.001666270       0.000035966       2.1584676  0.13395183E+02  0.17451136E+00          0          0
     2     _Null_     -0.093724464539    -284.255564723    -263.124994917     -21.130569805      -5.8696027  0.42804138E+01  0.22737545E+00          0          0
     3     Ztcd        0.024884736881     277.292523785     272.459894125       4.832629660       1.3711148  0.19779574E-04  0.18261173E+00          0          0
     4     Rad_a       0.017456962810     200.142269443     198.754025093       1.388244350       0.6984736  0.14990393E+01  0.17450967E+00          0          0
 Initial best average relative RMS error (%) =   69.9432080830336     
 Current best average relative RMS error (%) =   69.3784626874222     
 qmax =   13.3951834080603     
 qmed =   2.88972658510548     
 qmin =  1.977957368110705E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947182502308  0.48525432D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319548728317    0.302670964096  0.44354916D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001476618    0.000001398627  0.51493159D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.111908833267    0.105998088724  0.12574849D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 15
 qa,qb   =  -133.951636284866        133.951832102646     
 
 Summary Table - Ridge Regression Iteration   15
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95183408     0.00101279     0.00050640    -0.26117264    21.23193635      69.93004831         T & T         T & T
     2   108.24390671     0.00154598     0.00077299    -0.26203022    21.09773318      69.92392891         T & T         T & T
     3    85.24207695     0.00247885     0.00123942    -0.26352575    20.86248194      69.91335283         T & T         T & T
     4    64.94634482     0.00422549     0.00211274    -0.26630793    20.42059658      69.89398964         T & T         T & T
     5    47.35671030     0.00778072     0.00389036    -0.27188821    19.51649198      69.85628767         T & T         T & T
     6    32.47317340     0.01579645     0.00789823    -0.28395661    17.46642128      69.77910982         T & T         T & T
     7    20.29573411     0.03628257     0.01814128    -0.31020094    12.32338113      69.62266450         T & T         T & T
     8    10.82439245     0.10006430     0.05003215    -0.33908348     0.38109188      69.37846269         T & T         T & T
     9     4.05914841     0.38130710     0.19065355    -0.04661110    20.05610522      70.72677976         T & F         F & F
    10     3.00000000     0.52304288     0.26152144     0.19200806    26.00049916      71.55709767         T & F         F & F
    11     1.00000000     1.23189235     0.61594617     0.99510735    46.87738121      70.58658592         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.807434218546155     
 Total reduction in av_qd_rel_rmse (%) =   36.8326592149208     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.167471931067     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63412253E-03  0.25181790E-01  0.17022364E-02  0.51894004E-04       0.0304858  0.10000000E+01        1            0
     2         _Null_     0.12006074E-02  0.34649783E-01 -0.28425556E+03  0.78119374E+01       0.0274821  0.31954873E+00        1            0
     3         Ztcd       0.72648778E-03  0.26953437E-01  0.27729252E+03  0.52343724E+01       0.0188767  0.14766184E-05        0            0
     4         Rad_a      0.63410441E-03  0.25181430E-01  0.20014227E+03  0.20025235E+01       0.0100055  0.11190883E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25181790E-01  0.17022364E-02  0.51894004E-04       1.0000000     485.2543245     0.174511363     0.069758613     2.290606112
     2         _Null_     0.34649783E-01 -0.28425556E+03  0.78119374E+01       0.3195487       0.0044355     0.227375450     0.035003715     0.677059757
     3         Ztcd       0.26953437E-01  0.27729252E+03  0.52343724E+01       0.0000015       0.0051493     0.182611735     0.064143953     1.923527299
     4         Rad_a      0.25181430E-01  0.20014227E+03  0.20025235E+01       0.1119088       0.0125748     0.174509673     0.069759636     2.290684085
 Total relative Backus-Gilbert R(spread) =   7.18187725371871     

 *** Ridge Regression Iteration   16.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8242603677549     
 Optimum ridge parameter beta =   117.164612508948     
 Damping factor range qa,qb   =  -133.950001620480        133.950197558690     
 kbest, beta_best =           8   117.164612508948     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.017097240541       0.001737937       0.001702236       0.000035701       2.0972707  0.13395020E+02  0.17451984E+00          0          0
     2     _Null_     -0.093799552938    -305.517431498    -284.255564723     -21.261866776      -5.9060741  0.42801493E+01  0.22737324E+00          0          0
     3     Ztcd        0.024790709485     282.123610159     277.292523785       4.831086374       1.3521375  0.19791738E-04  0.18264044E+00          0          0
     4     Rad_a       0.017101597726     201.500583379     200.142269443       1.358313935       0.6786742  0.15012537E+01  0.17451814E+00          0          0
 Initial best average relative RMS error (%) =   69.3784626874222     
 Current best average relative RMS error (%) =   68.8276644262571     
 qmax =   13.3950199537864     
 qmed =   2.89070148429309     
 qmin =  1.979173840398442E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947170944345  0.47890753D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319532877450    0.302652257283  0.44116330D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001477545    0.000001399487  0.51314979D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112075510483    0.106154667103  0.12590313D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 16
 qa,qb   =  -133.950001620480        133.950197558690     
 
 Summary Table - Ridge Regression Iteration   16
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.95019954     0.00100971     0.00050485    -0.26815508    21.31574423      69.36563467         T & T         T & T
     2   108.24258586     0.00154130     0.00077065    -0.26902790    21.18142075      69.35960948         T & T         T & T
     3    85.24103679     0.00247140     0.00123570    -0.27055004    20.94596992      69.34919712         T & T         T & T
     4    64.94555231     0.00421302     0.00210651    -0.27338174    20.50374862      69.33013672         T & T         T & T
     5    47.35613243     0.00775861     0.00387930    -0.27906160    19.59911537      69.29303633         T & T         T & T
     6    32.47277714     0.01575545     0.00787772    -0.29134747    17.54863721      69.21714652         T & T         T & T
     7    20.29548646     0.03621033     0.01810517    -0.31809460    12.40938811      69.06364702         T & T         T & T
     8    10.82426037     0.09998858     0.04999429    -0.34813982     0.29281525      68.82766443         T & T         T & T
     9     4.05909887     0.38068433     0.19034217    -0.06008413    19.80050122      70.22955122         T & F         F & F
    10     3.00000000     0.52125652     0.26062826     0.17675823    25.70748193      71.09255567         T & F         F & F
    11     1.00000000     1.21820582     0.60910291     0.98232649    46.63412802      70.24436790         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.793903813704711     
 Total reduction in av_qd_rel_rmse (%) =   37.3341471424294     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.164612508948     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63420401E-03  0.25183407E-01  0.17379369E-02  0.52585115E-04       0.0302572  0.10000000E+01        1            0
     2         _Null_     0.12005886E-02  0.34649511E-01 -0.30551743E+03  0.78541238E+01       0.0257076  0.31953288E+00        1            0
     3         Ztcd       0.72680144E-03  0.26959255E-01  0.28212361E+03  0.52536815E+01       0.0186219  0.14775445E-05        0            0
     4         Rad_a      0.63418586E-03  0.25183047E-01  0.20150058E+03  0.20001923E+01       0.0099265  0.11207551E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25183407E-01  0.17379369E-02  0.52585115E-04       1.0000000     478.9075282     0.174519835     0.069756396     2.290310921
     2         _Null_     0.34649511E-01 -0.30551743E+03  0.78541238E+01       0.3195329       0.0044116     0.227373239     0.035008146     0.677158620
     3         Ztcd       0.26959255E-01  0.28212361E+03  0.52536815E+01       0.0000015       0.0051315     0.182640439     0.064127500     1.922429493
     4         Rad_a      0.25183047E-01  0.20150058E+03  0.20001923E+01       0.1120755       0.0125903     0.174518143     0.069757420     2.290388987
 Total relative Backus-Gilbert R(spread) =   7.18028802159632     

 *** Ridge Regression Iteration   17.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8241314599703     
 Optimum ridge parameter beta =   117.161821862719     
 Damping factor range qa,qb   =  -133.948406390123        133.948602325243     
 kbest, beta_best =           8   117.161821862719     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.016805683548       0.001773498       0.001737937       0.000035561       2.0461698  0.13394860E+02  0.17452809E+00          0          0
     2     _Null_     -0.093898058569    -326.928176133    -305.517431498     -21.410744635      -5.9474291  0.42798924E+01  0.22737111E+00          0          0
     3     Ztcd        0.024621093921     286.941119828     282.123610159       4.817509669       1.3303495  0.19791426E-04  0.18266842E+00          0          0
     4     Rad_a       0.016810043722     202.834801395     201.500583379       1.334218017       0.6621410  0.15034078E+01  0.17452640E+00          0          0
 Initial best average relative RMS error (%) =   68.8276644262571     
 Current best average relative RMS error (%) =   68.2883074537333     
 qmax =   13.3948604304385     
 qmed =   2.89165009113565     
 qmin =  1.979142620217084E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947159664341  0.47258564D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319517507945    0.302634095576  0.43855578D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001477539    0.000001399465  0.51107513D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112237657576    0.106306982076  0.12599173D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 17
 qa,qb   =  -133.948406390123        133.948602325243     
 
 Summary Table - Ridge Regression Iteration   17
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94860430     0.00100791     0.00050395    -0.27451569    21.35203751      68.81501985         T & T         T & T
     2   108.24129679     0.00153857     0.00076928    -0.27540427    21.21765544      68.80904962         T & T         T & T
     3    85.24002164     0.00246706     0.00123353    -0.27695394    20.98210908      68.79873340         T & T         T & T
     4    64.94477886     0.00420576     0.00210288    -0.27983707    20.53973309      68.77985298         T & T         T & T
     5    47.35556846     0.00774575     0.00387287    -0.28562097    19.63488411      68.74311887         T & T         T & T
     6    32.47239042     0.01573164     0.00786582    -0.29813730    17.58441937      68.66805574         T & T         T & T
     7    20.29524476     0.03616850     0.01808425    -0.32543560    12.44826985      68.51669923         T & T         T & T
     8    10.82413146     0.09994024     0.04997012    -0.35682560     0.26023449      68.28830745         T & T         T & T
     9     4.05905053     0.37998106     0.18999053    -0.07346706    19.62621956      69.74501691         T & F         F & F
    10     3.00000000     0.51935389     0.25967694     0.16152049    25.48746643      70.64150362         T & F         F & F
    11     1.00000000     1.20432544     0.60216272     0.96934493    46.42856439      69.91888398         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.783633989355659     
 Total reduction in av_qd_rel_rmse (%) =   37.8252180651409     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.161821862719     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63428337E-03  0.25184983E-01  0.17734980E-02  0.53291893E-04       0.0300490  0.10000000E+01        1            0
     2         _Null_     0.12005706E-02  0.34649251E-01 -0.32692818E+03  0.79007625E+01       0.0241667  0.31951751E+00        1            0
     3         Ztcd       0.72710718E-03  0.26964925E-01  0.28694112E+03  0.52761176E+01       0.0183875  0.14775388E-05        0            0
     4         Rad_a      0.63426521E-03  0.25184622E-01  0.20283480E+03  0.19989107E+01       0.0098549  0.11223766E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25184983E-01  0.17734980E-02  0.53291893E-04       1.0000000     472.5856354     0.174528091     0.069754241     2.290023517
     2         _Null_     0.34649251E-01 -0.32692818E+03  0.79007625E+01       0.3195175       0.0043856     0.227371111     0.035012453     0.677254612
     3         Ztcd       0.26964925E-01  0.28694112E+03  0.52761176E+01       0.0000015       0.0051108     0.182668423     0.064111468     1.921360071
     4         Rad_a      0.25184622E-01  0.20283480E+03  0.19989107E+01       0.1122377       0.0125992     0.174526398     0.069755267     2.290101626
 Total relative Backus-Gilbert R(spread) =   7.17873982715150     

 *** Ridge Regression Iteration   18.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8240060732284     
 Optimum ridge parameter beta =   117.159107473286     
 Damping factor range qa,qb   =  -133.946854809942        133.947050673284     
 kbest, beta_best =           8   117.159107473286     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.016547011940       0.001808979       0.001773498       0.000035481       2.0006004  0.13394705E+02  0.17453611E+00          0          0
     2     _Null_     -0.094005593120    -348.496100269    -326.928176133     -21.567924136      -5.9910900  0.42796433E+01  0.22736905E+00          0          0
     3     Ztcd        0.024409815376     291.738088930     286.941119828       4.796969102       1.3072858  0.19784176E-04  0.18269563E+00          0          0
     4     Rad_a       0.016551372873     204.147948619     202.834801395       1.313147224       0.6473974  0.15054981E+01  0.17453442E+00          0          0
 Initial best average relative RMS error (%) =   68.2883074537333     
 Current best average relative RMS error (%) =   67.7591647823136     
 qmax =   13.3947052651702     
 qmed =   2.89257066662851     
 qmin =  1.978417593630159E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947148692500  0.46636776D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319502607277    0.302616476732  0.43585833D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001477015    0.000001398952  0.50885913D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112395013368    0.106454789955  0.12604354D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 18
 qa,qb   =  -133.946854809942        133.947050673284     
 
 Summary Table - Ridge Regression Iteration   18
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94705265     0.00100675     0.00050338    -0.28055270    21.36394765      68.27577253         T & T         T & T
     2   108.24004293     0.00153681     0.00076841    -0.28145734    21.22953799      68.26983732         T & T         T & T
     3    85.23903423     0.00246428     0.00123214    -0.28303511    20.99394826      68.25958292         T & T         T & T
     4    64.94402655     0.00420111     0.00210055    -0.28597079    20.55150792      68.24082013         T & T         T & T
     5    47.35501989     0.00773750     0.00386875    -0.29186138    19.64659709      68.20433300         T & T         T & T
     6    32.47201426     0.01571637     0.00785818    -0.30461568    17.59634060      68.12986343         T & T         T & T
     7    20.29500966     0.03614148     0.01807074    -0.33249237    12.46284794      67.98024885         T & T         T & T
     8    10.82400607     0.09990316     0.04995158    -0.36532770     0.24905224      67.75916478         T & T         T & T
     9     4.05900351     0.37922737     0.18961369    -0.08684492    19.49441325      69.27145355         T & F         F & F
    10     3.00000000     0.51738289     0.25869144     0.14621456    25.30498828      70.20187196         T & F         F & F
    11     1.00000000     1.19035846     0.59517923     0.95610125    46.23949014      69.60640648         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.774865699780491     
 Total reduction in av_qd_rel_rmse (%) =   38.3069891242674     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.159107473286     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63436045E-03  0.25186513E-01  0.18089786E-02  0.54005691E-04       0.0298542  0.10000000E+01        1            0
     2         _Null_     0.12005531E-02  0.34648999E-01 -0.34849610E+03  0.79496012E+01       0.0228112  0.31950261E+00        1            0
     3         Ztcd       0.72740446E-03  0.26970437E-01  0.29173809E+03  0.53001774E+01       0.0181676  0.14770147E-05        0            0
     4         Rad_a      0.63434228E-03  0.25186153E-01  0.20414795E+03  0.19982105E+01       0.0097881  0.11239501E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25186513E-01  0.18089786E-02  0.54005691E-04       1.0000000     466.3677648     0.174536111     0.069752153     2.289744505
     2         _Null_     0.34648999E-01 -0.34849610E+03  0.79496012E+01       0.3195026       0.0043586     0.227369053     0.035016636     0.677347792
     3         Ztcd       0.26970437E-01  0.29173809E+03  0.53001774E+01       0.0000015       0.0050886     0.182695634     0.064095882     1.920320811
     4         Rad_a      0.25186153E-01  0.20414795E+03  0.19982105E+01       0.1123950       0.0126044     0.174534417     0.069753179     2.289822630
 Total relative Backus-Gilbert R(spread) =   7.17723573792033     

 *** Ridge Regression Iteration   19.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8238843380523     
 Optimum ridge parameter beta =   117.156472163534     
 Damping factor range qa,qb   =  -133.945348463289        133.945544214496     
 kbest, beta_best =           8   117.156472163534     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.016305280488       0.001844402       0.001808979       0.000035424       1.9582233  0.13394555E+02  0.17454389E+00          0          0
     2     _Null_     -0.094114830050    -370.224683540    -348.496100269     -21.728583271      -6.0357176  0.42794016E+01  0.22736706E+00          0          0
     3     Ztcd        0.024174731694     296.510230979     291.738088930       4.772142049       1.2837377  0.19772849E-04  0.18272206E+00          0          0
     4     Rad_a       0.016309640875     205.441606706     204.147948619       1.293658087       0.6336865  0.15075239E+01  0.17454220E+00          0          0
 Initial best average relative RMS error (%) =   67.7591647823136     
 Current best average relative RMS error (%) =   67.2396497921161     
 qmax =   13.3945546191781     
 qmed =   2.89346277475612     
 qmin =  1.977284925940919E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947138040219  0.46029115D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319488159858    0.302599389601  0.43313836D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001476186    0.000001398152  0.50658030D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112547522853    0.106598040227  0.12607381D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 19
 qa,qb   =  -133.945348463289        133.945544214496     
 
 Summary Table - Ridge Regression Iteration   19
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94554619     0.00100592     0.00050296    -0.28642155    21.36332845      67.74670252         T & T         T & T
     2   108.23882558     0.00153554     0.00076777    -0.28734244    21.22890717      67.74079228         T & T         T & T
     3    85.23807557     0.00246226     0.00123113    -0.28894864    20.99330082      67.73058240         T & T         T & T
     4    64.94329614     0.00419772     0.00209886    -0.29193755    20.55084230      67.71190585         T & T         T & T
     5    47.35448731     0.00773150     0.00386575    -0.29793636    19.64594738      67.67560572         T & T         T & T
     6    32.47164906     0.01570521     0.00785261    -0.31093295    17.59599227      67.60161302         T & T         T & T
     7    20.29478140     0.03612134     0.01806067    -0.33940333    12.46488770      67.45354299         T & T         T & T
     8    10.82388434     0.09986934     0.04993467    -0.37374009     0.24532946      67.23964979         T & T         T & T
     9     4.05895786     0.37844044     0.18922022    -0.10025591    19.38495089      68.80797069         T & F         F & F
    10     3.00000000     0.51537098     0.25768549     0.13080492    25.14172803      69.77256784         T & F         F & F
    11     1.00000000     1.17636887     0.58818444     0.94256246    46.05590590      69.30490808         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.766708078333800     
 Total reduction in av_qd_rel_rmse (%) =   38.7799944224190     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.156472163534     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63443519E-03  0.25187997E-01  0.18444024E-02  0.54721879E-04       0.0296692  0.10000000E+01        1            0
     2         _Null_     0.12005362E-02  0.34648754E-01 -0.37022468E+03  0.79994655E+01       0.0216071  0.31948816E+00        1            0
     3         Ztcd       0.72769313E-03  0.26975788E-01  0.29651023E+03  0.53250763E+01       0.0179592  0.14761856E-05        0            0
     4         Rad_a      0.63441702E-03  0.25187636E-01  0.20544161E+03  0.19978484E+01       0.0097247  0.11254752E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25187997E-01  0.18444024E-02  0.54721879E-04       1.0000000     460.2911547     0.174543890     0.069750131     2.289474072
     2         _Null_     0.34648754E-01 -0.37022468E+03  0.79994655E+01       0.3194882       0.0043314     0.227367057     0.035020697     0.677438243
     3         Ztcd       0.26975788E-01  0.29651023E+03  0.53250763E+01       0.0000015       0.0050658     0.182722059     0.064080748     1.919312146
     4         Rad_a      0.25187636E-01  0.20544161E+03  0.19978484E+01       0.1125475       0.0126074     0.174542195     0.069751157     2.289552199
 Total relative Backus-Gilbert R(spread) =   7.17577666011631     

 *** Ridge Regression Iteration   20.
 Return status code irr_lsstat =           0
 Optimum damping factor =   10.8237662755514     
 Optimum ridge parameter beta =   117.153916387763     
 Damping factor range qa,qb   =  -133.943887594549        133.944083208237     
 kbest, beta_best =           8   117.153916387763     
 Last realizable trial model  =          11
     j    pmxs_name    New_xp_best(j)  New pmxs_best(j)  Old pmxs_best(j)      pmxs(j)_step  percent change        s_svd(j)     resol_xp(j)  is_insens  is_freeze
     1     kiso        0.016072290400       0.001879775       0.001844402       0.000035373       1.9178322  0.13394409E+02  0.17455143E+00          0          0
     2     _Null_     -0.094222077549    -392.114883251    -370.224683540     -21.890199711      -6.0806110  0.42791673E+01  0.22736512E+00          0          0
     3     Ztcd        0.023925363739     301.254694210     296.510230979       4.744463231       1.2601154  0.19758958E-04  0.18274770E+00          0          0
     4     Rad_a       0.016076649430     206.716614774     205.441606706       1.275008069       0.6206182  0.15094860E+01  0.17454973E+00          0          0
 Initial best average relative RMS error (%) =   67.2396497921161     
 Current best average relative RMS error (%) =   66.7294979335018     
 qmax =   13.3944085184133     
 qmed =   2.89432662972323     
 qmin =  1.975895840059227E-005
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   1    1.000000000000    0.947127709335  0.45437190D+03
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   2    0.319474149061    0.302582818992  0.43043042D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   3    0.000001475165    0.000001397169  0.50427967D-02
 i,[qsv(i)/qmax],uqsvd(i),wjscale(i) =   4    0.112695233521    0.106736778378  0.12609057D-01
 Damping range setting =           1
 ntimes_qdset(1:iqdsetting) = 20
 qa,qb   =  -133.943887594549        133.944083208237     
 
 Summary Table - Ridge Regression Iteration   20
 -----------------------------------------------
     k        qdls(k)         rxp(k)      xp_rms(k)     sumXp/|Xp|  dtheta_pbm(k)       av_qd_rel_  is_xp_realiz(k) &  Search &
                                                                                        rmse_ls(k)  is_pmxs_realiz(k)   new_tmk
     1   133.94408518     0.00100523     0.00050261    -0.29220233    21.35629132      67.22724189         T & T         T & T
     2   108.23764497     0.00153449     0.00076725    -0.29313957    21.22186677      67.22135166         T & T         T & T
     3    85.23714584     0.00246059     0.00123030    -0.29477439    20.98625780      67.21117779         T & T         T & T
     4    64.94258777     0.00419492     0.00209746    -0.29781692    20.54380514      67.19257204         T & T         T & T
     5    47.35397079     0.00772652     0.00386326    -0.30392485    19.63896609      67.15642953         T & T         T & T
     6    32.47129487     0.01569589     0.00784794    -0.31716615    17.58935716      67.08285694         T & T         T & T
     7    20.29456004     0.03610404     0.01805202    -0.34623882    12.46047603      66.93623809         T & T         T & T
     8    10.82376628     0.09983486     0.04991743    -0.38210874     0.24416618      66.72949793         T & T         T & T
     9     4.05891359     0.37763054     0.18881527    -0.11371445    19.28737368      68.35410638         T & F         F & F
    10     3.00000000     0.51333470     0.25666735     0.11527881    24.98819382      69.35300245         T & F         F & F
    11     1.00000000     1.16239742     0.58119871     0.92871025    45.87202011      69.01324327         F & F         F & F
 
 Stage reduction in av_qd_rel_rmse (%) =  0.758706894208283     
 Total reduction in av_qd_rel_rmse (%) =   39.2444748253708     

 The scaled/raw parameter std deviations for the damped LS solution:
 Lagrange multiplier (betasq) =   117.153916387763     

 parameter  pmxs_name         var[xp(i)]    stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]  pmxs_std_ratio   qsv(i)/qsv(1)  ipmresflag(i)  iwjflag(i)
     1         kiso       0.63450760E-03  0.25189434E-01  0.18797750E-02  0.55437922E-04       0.0294918  0.10000000E+01        1            0
     2         _Null_     0.12005197E-02  0.34648516E-01 -0.39211488E+03  0.80497369E+01       0.0205290  0.31947415E+00        1            0
     3         Ztcd       0.72797324E-03  0.26980979E-01  0.30125469E+03  0.53503999E+01       0.0177604  0.14751647E-05        0            0
     4         Rad_a      0.63448942E-03  0.25189074E-01  0.20671661E+03  0.19976969E+01       0.0096639  0.11269523E+00        1            0

 parameter  pmxs_name       stdev[xp(i)]         pmxs(i)  stdev[pmxs(i)]   qsv(i)/qsv(1)      wjscale(i)      resol(i,i)   bg_rspread(i)  rel_bg_rspread
     1         kiso       0.25189434E-01  0.18797750E-02  0.55437922E-04       1.0000000     454.3719040     0.174551426     0.069748175     2.289212199
     2         _Null_     0.34648516E-01 -0.39211488E+03  0.80497369E+01       0.3194741       0.0043043     0.227365116     0.035024639     0.677526060
     3         Ztcd       0.26980979E-01  0.30125469E+03  0.53503999E+01       0.0000015       0.0050428     0.182747702     0.064066063     1.918333836
     4         Rad_a      0.25189074E-01  0.20671661E+03  0.19976969E+01       0.1126952       0.0126091     0.174549731     0.069749202     2.289290318
 Total relative Backus-Gilbert R(spread) =   7.17436241373262     
